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SUMMARY 

In recent years, simulation models have been used as a complementary tool for research and for 
quantifying soil carbon sequestration under widely varying conditions. This has improved the 
understanding and prediction of soil organic carbon (SOC) dynamics and crop yield responses to soil 
and climate conditions and crop management scenarios. The goal of the present study was to estimate 
the changes in SOC for different cropping systems in West Africa using a simulation model. A crop 
rotation experiment conducted in Farako-Ba, Burkina Faso was used to evaluate the performance of 
the cropping system model (CSM) of the Decision Support System for Agrotechnology Transfer 
(DSSAT) for simulating yield of different crops. Eight crop rotations that included cotton, sorghum, 
peanut, maize and fallow, and three different management scenarios, one without N (control), one 
with chemical fertilizer (N) and one with manure applications, were studied. The CSM was able to 
simulate the yield trends of various crops, with inconsistencies for a few years. The simulated SOC 
increased slightly across the years for the sorghum-fallow rotation with manure application. However, 
SOC decreased for all other rotations except for the continuous fallow (native grassland), in which the 
SOC remained stable. The model simulated SOC for the continuous fallow system with a high degree 
of accuracy normalized root mean square error (RMSE) = O001, while for the other crop rotations the 
simulated SOC values were generally within the standard deviation (s.d.) range of the observed data. 
The crop rotations that included a supplemental N-fertilizer or manure application showed an increase 
in the average simulated aboveground biomass for all crops. The incorporation of this biomass into the 
soil after harvest reduced the loss of SOC. In the present study, the observed SOC data were used for 
characterization of production systems with different SOC dynamics. Following careful evaluation of 
the CSM with observed soil organic matter (SOM) data similar to the study presented here, there 
are many opportunities for the application of the CSM for carbon sequestration and resource 
management in Sub-Saharan Africa. 



INTRODUCTION 

Several studies have shown that soils in West Africa 
are structurally degraded and that the nutrients are 
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depleted (Pichot et al. 1981; Bado et al. 1997; de 
Ridder et al. 2004; Schlecht et al. 2007). Furthermore, 
agricultural systems and farming techniques contrib- 
ute to the degradation of environmental resources 
(Bandre & Batta 1998). In recent years, interest in the 
sequestration and turnover of organic matter (OM) in 
soils has increased because soil organic matter (SOM) 
affects the stabilization of soil structure and the 
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cycling of organically held plant nutrients (Oades & 
Waters 1991; Tiessen et al. 1994; Lai et al. 2007). In 
addition, carbon (C) held in or released from SOM 
could change atmospheric carbon dioxide (C0 2 ) levels 
(Schlesinger 1991; Reeves 1997; Tiessen et al. 1998). 
The exchange of C between the land surface and the 
atmosphere is responsible for a large fraction of the 
global variation in atmospheric C0 2 , both within and 
between years (Tucker et al. 1986). Changes in 
vegetation productivity and C storage due to climate 
change and change in land use could have an impor- 
tant impact on the global C cycle (Mooney & Koch 
1994). 

Sorghum (Sorghum bicolor) is the main cereal crop 
in West Africa. Cotton (Gossypium spp.) is considered 
to be more profitable economically and plays an im- 
portant role in Burkina Faso's economy, but can only 
be grown in areas where there is sufficient rainfall. 
Peanut (Arachis hypogoea) and cowpea (Vigna un- 
guiculata) are the most commonly grown legumes 
(Bado 2002). In many cases, food crops are grown as 
dual purpose crops in which the seeds are used for 
human consumption, whereas the stover is used for 
feeding animals. Bado (2002) reported that rotations 
with peanut or cowpea increased the amount of avail- 
able nitrogen (N) for the next crop, resulting in a 
considerable increase in sorghum yield. Short-term 
use of chemical fertilizer, followed by maintenance 
fertilization using chemical and animal manure, are 
required to improve soil fertility and increase produc- 
tivity, but in some cases socio-economic constraints 
have largely prevented the improvement of most semi- 
arid agricultural systems to date (Tiessen et al. 1998; 
Bationo et al. 2007; Doraiswamy et al. 2007). 

Bado et al. (2006/?) emphasized the importance of 
identifying a cost-effective way of improving soil 
fertility and productivity in Burkina Faso and other 
countries in Sub-Saharan Africa. In a crop rotation 
study conducted in Burkina Faso, Bado et al. (2006ft) 
found that high yields were achieved by simultaneous 
application of organic and mineral fertilizers or when 
mineral fertilizers were associated with dolomite. 
Nitrogen fixing legumes, such as peanut and cowpea, 
are important components of the cropping system 
of resource-poor farmers and serve as food as well 
as cash crops (Bado et al. 2006a). Crop rotations of 
continuous sorghum and cotton-peanut-sorghum re- 
sulted in a decrease in soil organic carbon (SOC). Low 
quantities of manure, e.g. 5 tonnes(t)/ha every 2 years, 
did not significantly increase SOC after 6 years of 
cultivation. Ouedraogo et al. (2006) reported that com- 
bining recalcitrant organic amendments and N ferti- 
lizer is the best option in sustaining crop production in 
semi-arid West Africa due to a positive change in soil 
particulate OM and a reduction in the decline of soil 
C, especially for the more stable fraction that is re- 
sponsible for the soil structure maintenance. Hien 
et al. (2006) recommended selection of crops with a 



good rooting system, the application of organic 
manures, which are resistant to rapid biodegradation, 
reduction of tillage to avoid creating conditions for 
rapid biodegradation and regular liming in order to re- 
duce the loss of soil C. A long-term rotation study was 
conducted from 1960 to 1978 in a degraded ferrugi- 
nous tropical soil of Burkina Faso: soil C was 
measured in 1969 and 1978 and it was found that the 
soil C in the 0-200 mm layer decreased from 2-9 to 
2-5 g/kg for the control treatment without fertilizers, 
from 2-9 to 2-4 g/kg for the treatment with appli- 
cations of inorganic fertilizers and from 3T to 2-5 g/kg 
for the treatment with inorganic fertilizers and crop 
residues (sorghum) incorporated into the soil. The 
soil C increased from 3-1 to 3-5 g/kg for the treatment 
with inorganic fertilizers and 5 t/ha/yr of manure and 
increased from 5-3 to 6-6 g/kg for the treatment 
with inorganic fertilizers and 40 t/ha/yr of manure 
(Pichot et al. 1981). Lai (2000) observed annual rates 
of soil C increase under no-till management ranging 
from 363 kg/ha/yr to more than 1000 kg/ha/yr (for one 
severely depleted soil) over a 3-year experiment to 
restore soil C in western Nigeria. 

In other regions, such as India, it was found that the 
total system productivity was greater in groundnut- 
groundnut and groundnut-chickpea systems in com- 
parison with other systems, but the SOC declined 
(Ghosh et al. 2006). In a study conducted in Spain, 
Herencia et al. (2008) reported that in general, treat- 
ments with manure resulted in increased soil nutrient 
values, which had an effect on several later crop 
cycles, and that the use of organic composts resulted 
in an increase in SOC and the storage of nutrients, 
which can provide long-term fertility benefits. 

However, soils, climate and management practices 
vary over space and time, creating an almost infinite 
combination of factors that interact and influence 
how much C is stored in soils. Quantifying soil C 
sequestration under widely varying conditions is com- 
plicated. Thus, simulation models can be used to 
complement information gained from experiments to 
help understand and predict SOC and yield responses 
to soil, climate and crop management. Typically, 
long-term effects can be assessed by the models to 
predict environmental risks and to evaluate alterna- 
tive management practices for alleviating these risks 
(Muchow et al. 1991; Thornton et al. 1995; Chipanshi 
et al. 1997; Nijbroek et al. 2003). The Decision Sup- 
port System for Agrotechnology Transfer (DSSAT) is 
a comprehensive decision support system (DSS) for 
assessing management options (Tsuji et al. 1994; 
Jones et al. 2003; Hoogenboom et al. 2004). The 
cropping system simulation model (CSM) included in 
DSSAT v4.0 (Jones et al. 2003; Hoogenboom et al. 
2004), is process-oriented, dynamic and simulates 
growth, development and yield for more than 25 dif- 
ferent crops. The Century model (Parton et al. 1993) 
simulates the long-term dynamics of C for different 
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plant-soil systems. It was developed to analyse a wide 
range of crop rotations and tillage practices and to 
determine the impact of crop management on the 
productivity and sustainability of agroecosystems. 
The Century model has been used successfully to 
simulate SOC across a variety of land uses and climate 
types (Kelly et al. 1997; Diagana et al. 2007; Ogle 
et al. 2007). Gijsman et al. (2002) incorporated a 
Century-based module into CSM to allow for more 
flexibility in handling different agricultural systems for 
long-term simulations of crop rotations. This capa- 
bility is particularly important to enable CSM to be 
used for predicting yields in low-input cropping sys- 
tems where soils tend to be deficient in OM and 
nutrients. However, even under high-input manage- 
ment systems where large amounts of soluble nutrients 
are supplied, hence, minimizing the role of OM in 
supplying nutrients, its role in enhancing soil pro- 
ductivity via the modification of soil processes cannot 
be overlooked (Porter et al. 2010). 

There are a number of uncertainties associated with 
simulation models and their use, and if one does not 
adequately address these uncertainties, the simulated 
results will be meaningless (Jones et al. 2004). These 
uncertainties are due to the fact that models are sim- 
plifications of reality; there are uncertainties in model 
parameters and in the data inputs. Thus, research is 
needed to ensure that models can reproduce the re- 
sponses that are observed in field experiments. How- 
ever, soil C measurements are also uncertain and 
costly, and the errors could be much larger than the 
annual changes in SOC. By combining measurements 
with model predictions, a more accurate estimate of 
SOC can be obtained (Koo et al. 2003). The goal of 
the present study was to determine the yield variability 
and SOC dynamics for different rotations and 



fertilizer levels of cropping systems in a semi-arid 
region of West Africa. Specific objectives included the 
simulation of eight crop rotations and three fertilizer 
treatments, and the analysis of observed and simu- 
lated yield trends and SOC dynamics for the different 
crop rotations and management treatments. 

MATERIALS AND METHODS 

Crop management 

A crop rotation study was initiated in Farako-ba 
(4°20'W, 11°6'N), Burkina Faso in 1993 (Bado 2002). 
The rotations consisted of various cropping sequences 
that included sorghum, cotton, peanut and maize, 
which are important food and cash crops in the re- 
gion, and native grassland as a fallow. The experiment 
was initiated on a 6-year native grassland fallow and 
was an incomplete-factorial, split-plot design ar- 
ranged in randomized blocks with eight crop rotations 
as whole plot factors, nine input levels as split-plot 
factors and four replications. The input levels were 
distinguished by the type of soil amendment that was 
added (Bado 2002; Bostick et al. 2007). 

The eight rotations plus three management treat- 
ments were analysed from 1993 to 2004 (Table 1). The 
crop rotations were: (1) continuous fallow, (2) con- 
tinuous sorghum, (3) continuous peanut, (4) continu- 
ous cotton, (5) peanut-cotton-sorghum, (6) peanut- 
sorghum-cotton, (7) maize-sorghum-cotton and (8) 
sorghum-fallow. The three management treatments 
were: (1) no supplemental N fertilizer (control treat- 
ment), (2) the application of inorganic N and (3) the 
application of manure. The recommended application 
rates of inorganic N, P and K prior to planting were 
37, 10 and 1 1 kg/ha, respectively, for sorghum; 44, 
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Fig. 1. Total precipitation for the cropping season (May- 
October) in Farako-ba, Burkina Faso. 



15 and 17 kg/ha, respectively, for cotton; 14, 10 
and 1 1 kg/ha, respectively, for peanut and 60, 12 and 
7 kg/ha, respectively, for maize. For the manure treat- 
ments, air-dried cow manure (184gC/kg, 18gN/kg, 
3 ■ 1 g P/kg and 1 -6 g K/kg) at a rate of 3000 kg/ha was 
applied each year at planting. 



Climate 

The climate of Burkina Faso is similar to other 
countries in the region and is characterized by a dry 
season, which lasts, on average, from mid-November 
to mid-April, and a wet season, which lasts from May 
to mid-September. Most of the country (centre and 
south) lies in the Sudanese climatic zone, while the 
north is influenced by the climate of the Sahel. Rain- 
fall varies from 400 mm in the north to 1000 mm in the 
southwest. However, the most important aspects of 
the rainfall are its temporal and spatial variability 
(Bandre & Batta 1998). 

According to the FAO eco-climatic classification, 
Farako-ba is located in the Sudanian climatic zone, 
defined by one rainy season per year (FAO 2005). The 
average seasonal cumulative precipitation and s.d. 
from May to October from 1993 to 2004 at the study 
site were 927 and 171mm, respectively. The total 
precipitation for the cropping season (May-October) 
was variable for each individual year (Fig. 1). The 
maximum total seasonal rainfall was 1 1 74 mm for 
1998 and the minimum was 639 mm for 2002. The 
cropping season is also characterized by high temp- 
eratures. For example, the average maximum and 
minimum temperatures for June were 32 and 22-3 °C, 
respectively, while for July they were 30 and 21-5 °C, 
respectively. However, the highest temperatures 
usually occur during the months of March and April, 
with average maximum temperature of 38 °C. 



Soil characteristics and measurements 

The soil at the experimental site was a weathered 
ferruginous tropical soil (Bado 2002), which is an 
Alfisol or Ultisol in the USDA soil taxonomy 
classification (Soil Survey Staff 1999). This soil covers 
0-40 of the land area in the sub-humid zone of West 
Africa (6-1 2°N; Adeoye & Mohamed-Saleem 1990). 
These ferruginous tropical soils are characterized by a 
sandy surface horizon that is low in OM and base 
exchange capacity (ILCA 1979). The soil texture at 
the experimental site was loamy sandy, with sand, silt 
and clay contents of 74, 19 and 7% by volume, 
respectively (Bado 2002; Bostick et al. 2007). The soil 
profile inputs were complemented with data from a 
ferruginous tropical soil from northern Ghana that 
was obtained from DSSAT V4.0 database (Gijsman 
et al. 2007). The next soil horizon used as input for the 
model was similar to the upper horizon in terms of 
texture, and was characterized by sand, silt and clay 
contents of 66, 21 and 13% by volume, respectively 
(0-2-0-40 m depth). The next soil horizon, 0-4-1 -2 m 
depth, contained 55, 22 and 23% by volume of sand, 
silt and clay, respectively. The bulk density in the top 
0-2 m of the soil was estimated as 1 -46 g/cm 3 from soil 
texture components (Rawls 1983). This top layer of 
the soil had low values for total OM (1 1 g/kg), a weak 
cation exchange capacity (CEC), was slightly acidic 
(pH water = 6-5) and was N and P poor (Bado 2002). 

Three composite samples were collected from the 
top 0-2 m of soil before planting in 1993 to character- 
ize the initial SOC of the native grassland. The 
average SOC was 5-5 g/kg. SOC was also measured 
after harvest in 1998 and 2003 for some crop rota- 
tions only (continuous fallow, continuous sorghum, 
sorghum-fallow, continuous cotton and maize- 
sorghum-cotton) and input levels (control, N-fertili- 
zer and manure applications). The Walkley-Black 
method (Walkley & Black 1934) was used for all SOC 
determinations. For the 1998 and 2003 SOC measure- 
ments, four composite soil samples were taken and 
analysed for each treatment that was measured. 



Simulation of crop rotations 

The crop rotation or sequence tool of CSM-DSSAT 
V4.0 (Thornton et al. 1995) was used to simulate the 
eight crop rotations and three fertilizer treatments. 
The CSM-DSSAT sequence analysis tool allows the 
user to conduct simulations of crop rotations or crop 
sequences and to analyse the results (Bowen et al. 
1998). The main aspect of the sequence analysis is the 
consideration of experiments that are conducted 
across multiple cropping seasons. Therefore, the 
carry-over of the soil water and nutrient status is 
affected from one cropping season or crop to the 
subsequent one (Thornton et al. 1994). The Century 
soil OM module is integrated into the CSM-DSSAT 
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V4.0 in order to better model the dynamics of soil 
organic nutrient processes (Gijsman et al. 2002; Porter 
et al. 2010). The current version of CSM does not 
consider the P and K balance in the soil (Jones et al. 
2003); thus, these elements were not defined as inputs. 

Daily rainfall data were obtained from a conven- 
tional weather station located at the research station 
of the Environmental and Agricultural Research 
Institute (INERA), located adjacent to the experimen- 
tal site, in Farako-ba. The daily maximum and mini- 
mum temperature and solar radiation were obtained 
from the nearest weather station (10 km), located 
in Bobo Dioulasso. Using weather data from these 
two stations was the best way to obtain accurate 
daily temperatures and also rainfall at the experimen- 
tal site. 

The soil parameters needed to characterize the pro- 
file were obtained from Bado (2002). The parameters 
related to the different SOC decomposition rates that 
were required to run the Century component of the 
CSM were set according to previous reports on similar 
soils from Burkina Faso (Hien et al. 2006; Barthes 
et al. 2008). Therefore, the SOC corresponding with 
the fraction of stable organic C (SOM3, also called 
passive) was set equal to 0-26. Hien et al. (2006) 
reported 2-6 g/kg of SOC for a soil in the Savannah 
region of the Southwest of Burkina Faso after 13 years 
of continuous cropping, with removal of the crop 
residues. The intermediate SOC (SOM2) was set equal 
to 0-72. The model was set up to initiate the simu- 
lations 6 years prior to the beginning of the exper- 
iment. During this period, the growth of a native 
grassland fallow was simulated, mimicking the same 
crop that was grown prior to the start of the rotation 
experiment. Thereafter, the eight rotations and the 
three management treatments were simulated from 
1993 to 2004 (Table 1). 

For model evaluation, the simulated dates of 
flowering and maturity as well as yield were compared 
with the observed values. The phenology was eval- 
uated using growth and development data for 2004 
(data not shown). The cultivar coefficients are specific 
for each crop, and were obtained sequentially, starting 
with the phenological development parameters related 
to flowering and maturity dates, followed by the crop 
growth parameters related with kernel filling rate and 
kernels number per plant (Hunt & Boote 1998). An 
iterative procedure (Hunt et al. 1993) was used to 
select the most appropriate value for each coefficient. 
A detailed description of the cultivar coefficients 
used by the CSM-DSSAT V4.0 can be found in 
Hoogenboom et al. (1994) and Hunt & Boote (1998). 
The combination of cultivar coefficients that resulted 
in an average simulated yield that was the most similar 
to the average observed yield was selected. The sel- 
ected varieties used in the present study were: Sariaso- 
02 for sorghum, RMP12 for peanut, FK37 for cotton 
and SR22 for maize. 



Statistical analysis 

An analysis of yield variability, total biomass incor- 
porated into the soil and SOC was conducted for the 
different crop rotations and fertilizer treatment com- 
binations. The observed and simulated yield vari- 
ations were examined by plotting the normalized yield 
for each crop in each rotation and treatment according 
to Meinke & Hammer (1995). 

Y n = (Yi-Y)/S.D. (1) 

where Y„ is the normalized yield, Yi is the yield for 
each individual growing season, Y is the average yield 
for each crop across all years of the rotation and s.d. is 
the standard deviation of the yield for each crop across 
all years of the rotation. 

For the simulated yield, an analysis of the trends 
was conducted for the different rotations. The results 
were examined by plotting the progression with time 
of the cumulative deviations from the mean for the 
different crops of the different rotations and treat- 
ments. The plots start and end with zero, and the year 
to year fluctuations correspond to periods above or 
below the average yield for different crops. In the pres- 
ent study, the plots show long-term yield trends for 
different crops without removing variability (Russell 
1981; Hammer et al. 1987; Meinke & Hammer 1995; 
Garcia et al. 2006). For the simulated SOC dynamic, 
an analysis of the different rotations and treatments 
was conducted, which was presented in a graphical 
format depicting also the observed average and s.d. of 
the SOC values for the top soil layer (0-0-2 m) for 
1993, 1998 and 2003. 

RESULTS 

Yield 

The crop rotation or sequence tool of DSSAT V4.0 
was able to simulate the observed yields for different 
crop rotations (Table 2). The average observed yields 
varied between 720 kg/ha for the continuous sorghum 
rotation and 2005 kg/ha for maize crop in the cotton- 
maize-sorghum rotation. The simulated average 
yields varied between 801 kg/ha for the continuous 
cotton rotation and 2672 kg/ha for the maize crop in 
the cotton-maize-sorghum rotation. The Pearson co- 
efficient between observed and simulated yield values 
was 0-84 and the root mean square error (RMSE, 
expressed as a proportion) was 0-26 indicating a fairly 
good simulation. The variation in the simulated yield 
was similar to the variation of the observed yield for 
all crop rotations and the three treatments (Figs 2-4). 
There were a few anomalies in which the simulated 
variation in yield differed from the observed, e.g. in 
1993 the models in general under-predicted the yield. 
The differences between simulated and observed 
values were partially due to the fact that the models 



584 C. M. TOJO SOLER E T A L 



Table 2. Observed and simulated average yields for different crop rotations 
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Fig. 2. Simulated and observed yield expressed as s.d. from the mean for the different rotations without N fertilizer (control). 



do not consider all biotic and abiotic factors that 
impact crop growth and development. Differences 
between observed and simulated values were probably 
due to pest and diseases, but there were no detailed 
experimental records to support this on particular 
years and crops. In general, the observed yields were 
higher than the average only at the beginning of the 
cropping system rotations. This can be explained by 
the better soil properties for the first year of the 
rotations, which were favoured by the 6 previous years 
with native grassland. However, for the remainder of 
the years for the different rotations, the variation in 
the simulated yield was, in general, similar to the 



observed yield. Lack of agreement between observed 
and simulated accumulated deviations occurred in 
the control treatments for the continuous sorghum 
and sorghum-fallow rotations in 2001, in which the 
observed yields were below the average, but the 
simulated yield values were above average (Fig. 2). 
For the peanut-cotton-sorghum rotation without N 
fertilizer (control), the simulated peanut yield for 1996 
and 2002 were below the average simulated, while the 
observed yields were above the average observed. In 
1999, the opposite occurred; the simulated yield was 
above the average and the observed yield was below 
average. These examples indicate that the model was 
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unable to simulate yield trends accurately for a few 
specific years for the control rotation (without N). 
Similar to these specific circumstances, there were also 
a few predictions in simulated yield trends that were 
different from observed for specific cases for crop 
rotations with N fertilizer and with manure appli- 
cations. This included the peanut-cotton-sorghum 
rotation with N fertilizer and also the rotation with 
manure applications for 1999 where the observed 
peanut yield was below the average, while the 
simulated yield was above the average (Figs 3 and 4). 

The progression across years of the cumulative 
deviations from the mean for the different crops of the 



various rotations and treatments were plotted for 
simulated and observed yield (Figs 5-7). The plots 
start and end with zero, and the year to year fluc- 
tuations correspond to periods above or below the 
average yield. Thus, the yield trends expressed as the 
cumulative deviation from the mean yield for each 
crop of a rotation and treatment had, for most of the 
years, crops and treatments, a good agreement be- 
tween the simulated and observed values (Figs 5-7). 
There was a decrease in observed yield for the final 
years of the crop rotations (2002^1), with a below- 
average yield for the different crops. In general, the 
decrease in yield for the final years of the rotation was 
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also simulated. The yield decrease can be explained by 
the reduction in soil fertility after years of continuous 
cropping systems. 

Biomass 

For the crop rotations involving N fertilizer and 
manure applications, the simulated average biomass 
and s.d. bars varied among rotations and crops. 
However, this variation was similar between N fer- 
tilized treatments and manure applied treatments, as 
represented in the average and s.d. bars of Fig. 8. 
When comparing average biomass values it was found 
that peanut was the crop with the largest simulated 
biomass values. The growing conditions were better 
for peanut because it is a legume and can fix N from the 
rhizobia that are formed in the root nodules. Sorghum 



was affected by N stress even when N was applied, 
because the amount was insufficient to supply the N 
needed by the crop. In the maize-sorghum-cotton 
rotation, the sorghum biomass with N fertilizers was 
lower than sorghum in the same sequence with 
manure. The manure applied increased the N available 
in the soil and this explains the high average values of 
biomass in comparison to the N fertilizer treatment. 
The aboveground biomass for the crop rotations 
without fertilizer (control) was removed from the 
experimental plot at harvest and not incorporated 
into the soil, causing a reduction in SOC across years. 

soc 

Both the observed and simulated SOC content for all 
crop rotations decreased across the years of the 
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experiment, except for the continuous fallow (native 
grassland) in which the SOC remained stable around 
5-0-5-5 g/kg (Fig. 9) and the sorghum-fallow rotation 
with manure application. The crop rotation or 
sequence tool of CSM-DSSAT V4.0 (and the inte- 
grated Century SOM module) was able to simulate the 
SOC for the continuous fallow system accurately, with 
a value of 0-01 for the normalized RMSE (expressed 
as a proportion). In general, the simulated SOC values 
were within the range of the s.d. of the observed data 
for the other crop rotations (Fig. 9). The observed 
SOC rate of decline in the rotations sorghum-fallow, 
continuous sorghum and continuous cotton appeared 
to stabilize after 1999, but with an important 
variability in the SOC values. The model predicted a 



continued reduction in SOC rate of decline that 
eventually could impact the sequestration of SOC 
estimations. The simulated biomass that was incor- 
porated into the soil for the N fertilized treatments 
and manure treatments resulted in a less pronounced 
decrease of simulated SOC across the years (Fig. 8). 

The crop rotation that had the best performance 
with respect to simulated SOC was the sorghum- 
fallow rotation with manure application because this 
rotation had an average increase of 3 1 kg/ha/yr for the 
top 0-2 m of the soil profile (Fig. 10). The same crop 
rotation showed smaller decreases in SOC for the N 
fertilizer application and for the treatment that did not 
have any fertilizer application in comparison with the 
other crop rotations. The decrease in simulated SOC 
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was 173kg/ha/yr for the N fertilizer application 
treatments and 232 kg/ha/yr for treatments that did 
not have any fertilizer. Decreases in SOC of more than 
400 kg/ha/yr were simulated for each of the other 
control crop rotations (treatments without fertilizer 
application) (Fig. 10). When N fertilizer was applied, 
the model simulated a decrease in SOC that ranged 
from 250 to 380 kg/ha/yr. When manure was applied, 



the model simulated a decrease in SOC decrease that 
ranged from 160 to 325 kg/ha/yr for the top 0-2 m of 
the soil, for continuous peanut and continuous cotton 
rotations, respectively. 

In cropping systems with limited amounts of 
external inputs, the SOM pool may contribute signifi- 
cantly to plant nutrition (Vanlauwe et al. 1998). The 
Century model includes three SOM pools, active 
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Table3. Simulated average SOC for the three components SO Ml (microbial), S0M2 (intermediate) and SO M3 
( passive) for all treatments at a depth of 0-0-2 m for 2004 after 11 years of crop rotation simulation 



SOM 1 (kg/ha) SOM 2 (kg/ha) SOM 3 (kg/ha) Total SOM (kg/ha) 



Control (without N) 154 6018 2836 9008 

Nitrogen fertilizer 341 6919 2867 10127 

Manure application 268 8185 2863 11316 

SOM 1 SOM 2 SOM 3 

(proportion) (proportion) (proportion) 

Control (without N) 0-02 0-67 0-32 

Nitrogen fertilizer 0-03 0-68 0-28 

Manure application 0-02 0-72 0-26 



(microbial SOM1), slow (intermediate SOM2) and 
passive (SOM3), which have different potential 
decomposition rates (Parton et al. 1993; Batjes 2001). 
At the end of crop rotation in 2004, the simulated 
average values for the different components of the 
SOC was high for the SOM2 fraction (intermediate) 
for the treatments with manure applied, accounting 
for 0-72 of the total (Table 3). On the other hand, the 
control crop rotations (without N applied) had the 
lowest values of SOM2 in year 2004. When analysing 
the SOC dynamic for each rotation, it was found that 
the continuous fallow (native grassland) had a high 
simulated SOM2 in year 2004 of 0-75. For the crop 
rotations without N fertilizer (control), the SOC 
decreased to near 0-003 after 12 years. 

DISCUSSION 

The removal of aboveground biomass at harvest for 
the crop rotations without fertilizer (control) caused 
an average reduction in SOC of 46% after 12 years, as 
there was no incorporation of biomass and residue 
into the soil. The removal of crop residues has been 
reported to be adverse in terms of soil quality due to 
the depletion of the soil C pool, which negatively 
affects crop production (Lai 2007). In a study 
conducted in Ghana, Adiku et al. (2009) found that 
over a 4-year period, maize rotated with bare fallow 
produced an average maize biomass and yield of 4-0 
and l-0t/ha/yr, respectively; SOC declined over time, 
with the maize bare fallow rotation losing about 55% 
of the initial SOC in the 4 years. The decline in SOC 
was 19% for a fertilized maize-grass rotation, and 
other treatments lost between 33 and 44% of the initial 
SOC. Furthermore, results from Zimbabwe reported 
by Mtambanengwe & Mapfumo (2005) highlight that 
the cumulative effects of applying substantial amounts 
of OM to specific areas on a regular basis resulted in 
consistently higher levels of SOC in designated rich 
fields than corresponding poor fields (or field sec- 
tions). The key challenge is enabling farmers to 
generate sufficient quantities of biomass during each 



cropping cycle in order to maintain the productivity of 
rich fields, while rehabilitating the non-productive 
ones (Mtambanengwe & Mapfumo 2005). Crop 
residues also have been reported to be a major factor 
in maintaining or reducing the loss of SOC on two 
tropical soils of Brazil (Leite et al. 2009). Nziguheba 
et al. (2005) found in western Kenya that organic 
residues can contribute to the improvement of soil 
nutrient cycling. However, the magnitude of their 
contribution depends on the biochemical properties of 
the residues. Regular applications of manure on rice- 
wheat and maize rotation also have been associated 
with an increase in SOC and grain yield across sites in 
China (Zhang et al. 2009). 

Zingore et al. (2005) reported that for soils on 
smallholdings in Zimbabwe, the organic C declined 
rapidly under cultivation to attain a new equilibrium 
within 10 years after woodland clearance. Mando 
et al. (2005) found that 10 years of continuous 
sorghum cultivation without organic inputs in a 
sandy loam Lixisol soil at Saria, Burkina Faso caused 
significant losses of C. However, addition of manure 
was effective in maintaining similar levels of C com- 
pared to fallow plots. Without manure, SOM was 
mainly stored as fine organic matter (FOM) with a 
size-fraction < 0-053 mm. In plots with manure and in 
fallow plots, the addition of manure doubled the par- 
ticulate organic matter (POM) (size-fraction 0-053- 
2 mm) concentration. POM was greatly affected by 
long-term soil management options. Similar results 
were also found by Nziguheba et al. (2005) in a soil 
amended with organic residues in western Kenya and 
by Ouattara et al. (2006) in western soils of Burkina 
Faso. Willson et al. (2001) reported that the POM 
amount combined with information about recently 
incorporated crop residues, were found to be good 
predictors of N mineralization potential in systems 
using conventional, legume-based organic and man- 
ure-based organic fertility management. In a study 
conducted in West Africa, Vanlauwe et al. (1998) 
found that the average proportion of soil N belonging 
to the POM pool ranged from 0-09 to 0-29 depending 
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on the annual N inputs from maize stover and 
prunings. Bationo & Buerkert (2001) recommended 
targeted applications of mineral fertilizers and effec- 
tive recycling of organic amendments as crop residues 
and manure to maintain food production for a rapidly 
growing population in the agro-pastoral farming 
systems in the Sudano-Sahelian zone of West Africa. 
Similarly, Vagen et al. (2005) and Farage et al. (2007) 
stated that for croplands the critical issue in increasing 
SOC is the addition of crop residues or manure and 
the reduction in soil disturbance through conservation 
tillage or no-till systems. De Costa & Sangakkara 
(2006) highlighted the benefits of agronomic measures 
such as mulching, the application of green manure, 
the use of crop residues, fallow rotation and agrofor- 
estry as very economical and efficient methods of 
increasing productivity and maintaining sustainability 
of small farming units in Asia. Zougmore (2003) 
recommended combining soil and water conservation 
techniques with soil fertility management practices to 
improve the efficiency of these measures, and thus 
sustain soil productivity in Burkina Faso and other 
countries in Sub-Saharan Africa. 

In the present study, the model simulated the SOC 
for the continuous fallow system accurately, with a 
normalized RMSE (expressed as a proportion) of 
0-001. In general, the simulated SOC values were 
within the range of the s.d. of the observed data for the 
other crop rotations. The simulated results for SOC 
from the present study are similar to previous studies 
that were conducted with the standalone version of the 
Century model. For instance, the Century model was 
able to accurately simulate the level of total soil C for 
five soil types in Brazil (Cerri et al. 2004): modelling 
also provided a flexible and powerful way to assess 
how different scenarios for pasture management and 
land use change can affect soil C dynamics (Cerri et al. 
2004). The present study highlights that the CSM is a 
useful tool to simulate crop growth and the develop- 
ment for various crop rotations and management 
practices. 

The understanding of the SOC dynamic for differ- 
ent crop rotations is important in relation to soil 
fertility aspects, but it also complements the research 
focus in recent years on the global carbon cycle, as it 
has become clear that increased levels of CO? and 
other greenhouse gasses in the atmosphere are causing 
changes to our climate at an increasingly rapid rate 



(McCarthy et al. 2001). Soil C sequestration has been 
recognized as an effective, low-cost technology to miti- 
gate climate change. Simulation models, alone or in 
combination with soil sampling and other techniques, 
have been used to help with monitoring changes in soil 
carbon levels as affected by climate, soil and crop 
management options (Apezteguia et al. 2009). How- 
ever, as stated by Batjes (2004), mitigation of climate 
change by increased C sequestration in soils appears 
particularly useful when addressed in combination 
with other pressing regional challenges that affect the 
livelihood of the people, such as combating land 
degradation and ensuring food security, while at the 
same time curtailing global anthropogenic emissions. 

In summary, the CSM simulated the yield trends of 
various crops grown in Burkina Faso fairly well, with 
some inconsistencies for a few years. The simulated 
SOC increased slightly across the years for the 
sorghum-fallow rotation with manure application. 
However, the most accurate simulation of SOC was 
obtained for the continuous fallow (native grass) with 
no supplemental fertilizer applied, in which the SOC 
remained stable around 5-0-5-5 g/kg. The crop 
rotations that included N-fertilizer or manure appli- 
cations showed an increase in the total amount of 
biomass that was simulated, which was incorporated 
back into the soil in order to reduce the loss of SOC. 
The present study showed that DSSAT and its 
associated CSM can be used for the identification of 
cropping systems that can enhance the agricultural 
production and sustainability in Burkina Faso and 
other countries in West Africa that have similarly 
degraded soils. However, it is very important to define 
the initial values for SOC accurately as this will affect 
the outcome of the simulations. Further research is 
needed in order to extend the adoption of CSMs as a 
tool for decision makers in Burkina Faso and other 
countries in the region aiming to recommend suitable 
rotations to different soils or environmental con- 
ditions. 
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